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ABSTRACT 


Design and analysis of control systems for flexible structures require accurate math 
models of flexible structures and software design with the analysis tools capable of 
Handling these models while maintaining numerical accuracy. Since aeroelastlc models 
; ‘ c of flexible structures tend to be large (e.g., 100 states), the availability of tools to handle 

sHch large models is crucial. Initial model development Is based on aerodynamic ^ 

, - mathematical models, wind tunnel data, mathematical structural models, and ground 

shake test results. Eventually, flight test data are used to update and refine the model. ~ 

i ^ This paper describes Boeing software tools used for the development of control laws of 

flexible structures. 

The Boeing Company has developed a software tool called Modem Control Software 
Package (MPAC). MPAC provides the environment necessary for linear model 
development, analysis, and controller design for large models of flexible structures. There 
are two features of MPAC which are particularly appropriate for use with large models: (1) 
numerical accuracy and (2) label-driven nature. With the first feature MPAC uses double 
precision arithmetic for all numerical operations and relies on EISPAC and LINPACK for 
the numerical foundation. With the second feature, all MPAC model inputs, outputs, and 
slates are referenced by user-defined labels. This feature allows model modification while 
maintaining the same state. Input, and output names. In addition, there Is no need for the 
user to keep track of a model variable's matrix row and column locations. 

There Is a wide range of model manipulation, analysis, and design features within the 
numerically robust and flexible environment provided by MPAC. Models can be built or 
modified using either state space or transfer function representations. Existing models 
can be combined via parallel, series, and feedback connections; and loops of a closed -loop 
model may be broken for analysis. Analysis tools available Include: 

eigenvalue/eigenvector, controllability matrix, observability matrix, transfer function 
generation, frequency response and singular value plots, covariance response to white 
noise or atmospheric turbulence models, model simulation using step, sinusoidal, 
random, or user-defined inputs. Control system design tools include: root locus LQG lull 
state feedback gain matrix computation, LQG full-order estimator design, and robust low 
order controller (SANDY) design as developed by Dr. Uy-Lot Ly at Stanford. 
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The existing Boeing Company structural analysis and design software package. ATLAS, 
has been extended in order to form a state-space model for Input to MPAC. The new 
capability, a module named DYFORM. is an outgrowth of earlier work under a NASA 
contract for Integrated Application of Active Controls. The structural and theoretical 
aerodynamic mathematical model originates within ATLAS In exactly the same fashion 
as for conventional flutter and dynamic loads analyses. The DYFORM module is then used 
to construct the state-variable model as required by MPAC. Its capabilities Include (1) 
control surfaces and/or gust vector as inputs, (2) sensors and/or loads quantities as 
outputs. (3) formulation in body-fixed or inertial axes. (4) modification of the theoretical 
aerodynamics using wind tunnel/flighl test data from rigid or flexible-model tests, and (5) 
use of S-plane rational airloads expressions to formulate the state model including 
augmented states to represent unsteady aerodynamic effects. 

MPAC has been used for yaw damper design (including active flexible mode suppression) of 
the Boeing 767 and 747 airplanes. The flexible structural models of these planes, as large 
as 100 states, have been handled by MPAC without loss of numerical accuracy. 

The Boeing Company plans for the development of a system identification and parameter 
estimation (SIPE) software tool. The system identification algorithms employ a multiple 
stepwise regression technique to determine the structure of the system. The parameter 
estimation algorithms update the current model using maximum likelihood estimation. 
The SIPE routines will be compatible with MPAC and RFJDATA (a data correction and 
reformatting program also developed by Boeing). The SIPE routines will be flexible, 
allowing the user to select gradient methods. Integration algorithms, and Riccati solution 
algorithms. The MPAC compatible model structure slated for the SIPE package will be 
applicable to any dynamic system. Aerodynamic, aeroelastic. ground effects, and sensor 
noise modeling will all be possible. 
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INTRODUCTION TO MPAC: 
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o MPAC was originally developed as a batch process tool. 
An interactive interface is currently being developed 
for MPAC to improve its ease of use and efficiency. 
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Provission for customized, user defined commands 
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COMMAND: 

♦COMPUTE THE OBSERVABILITY MATRIX OF THE Ug) 
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This is an example MPAC command 
command file is given on the 


file. The output file generated using this 
following pages. 


‘MPAC READ PLANT 
LAT2 . MDL 

‘COMPUTE EIGENVALUES OF PLANT 

*DEFINE PLANT 
DELETE STATE PS I 

CREATE STATE BETA_INT . dot : 1. BETA 

CREATE OUTPUT PHI_CRIT: 1. PHI. dot 5.0 PHI 

CREATE OUTPUT BETA_CRIT: 1. BETA. dot 3.2 BETA 4. BETA INT 


‘DESIGN GAIN MATRIX FOR PLANT 
.001 
2 , 2 

'AIL' 1. 

'RUD' 2. 

' PHI_CRIT' 4. 

' BETA_CRIT' 1. 

' WLOCUS' 'RHO' 1., 1., 1 

$ FIRST CUT LATERAL GAIN LOCUS 
$ AIL=1 . RUD=2 . 

$ PHI_CRIT=4 . BETA_CRIT=1. 

‘CREATE CONTROLLER FOR PLANT 
'NODIRECT' 

‘FORM PLANT + CONTROLLER 
‘PRINT CLOSED-LOOP SYSTEM 

‘COMPUTE EIGENVALUES OF CLOSED-LOOP SYSTEM 

‘MPAC WRITE CLOSED-LOOP SYSTEM 
CLOSED LOOP. MDL 
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MPAC output file example. Output file generated using command file on previous 

MPAC RELEASE VERSION 4.00 05 MAY 1987 (CONFIGURATION CONTROL I 

******************************************************* 

* MPAC INPUT/OUTPUT FILE DESCRIPTION === * 

* * 7 ****** *********************** *********************** 

COMMAND FILE example.cmd 

MODEL FILE 

INPUT BINARY FILE 

OUTPUT BINARY FILE 

MPAC OUTPUT FILE example. out 

MPAC GGP PLOT FILE example. ggp 

MPAC USER DATA FILE NO . 1 

MPAC USER DATA FILE NO. 2 

MPAC USER-DEFINED UBIN FILE — Dm (PST) 

TIME OF MPAC JOB EXECUTION Tuesday, July 5, 1988 3.50.36 pm ( b 

******************************************************* 

************************************************************************* 

************************************************************************* 

************************************************************************* 

*** ****** MODERN CONTROL THEORY ANALYSIS/SYNTHESIS SOFTWARE PACKAGE ****** 

****** APOLLO- VERS ION: MPAC 4.00 ON APOLLO FORTRAN 8.40 ****** 

I!*in******************^************************************************ 

***************************************************************:::::::::: 

************************************************************* 

************* 

* 07/05/88 * 

* 15:50:54 * 

************* 

********************************** 

************* TASK 1 ************* 

********************************** 

****** ****** 

****** *MPAC READ PLANT ****** 

****** ****** 

************************** ******** 

************* TASK 1 ************* 

********************************** 


*** MODEL READ FROM FILE: LAT2.MDL *** 

ELAPSED TIME (SEC): 0.24 
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******** 

****** 

****** 

****** 

******** 

******** 


****** 

♦COMPUTE EIGENVALUES OF PLANT ****** 

****** 


SAMPLING TIME : DELTA = 


0.0000 


******** EIGENVALUES OF PLANT *********************** 


COUNT 

REAL PART 

IMAG PART 

DAMPING 

FREQ (RAD/S) 

FREQ (HZ) 

1 

2 

3 

4 

5 

0.0000 
-1 . 2989E-02 
-0 . 1403 
-0 .1403 
-1.946 

0.0000 
0.0000 
1.676 
“1 . 676 
0.0000 

0.0000 

1.000 

8 . 3426E-02 
8 . 342 6E-02 
1.000 

0.0000 
1 .2989E-02 
1.682 
1.682 
1.946 

0.0000 

2 . 0673E-03 

0.2677 

0.2677 

0.3097 

ELAPSED 

TIME (SEC) : 

0.14 





*********** TASK 3 ************ 

****** * * * * * * 
****** *DEFINE PLANT ****** 
****** * * * * * * 

*********** TASK 3 ************ 




DELETED 

STATE 

PSI 



CREATED 

STATE 

BETA INT. : 

1.000 

BETA 

CREATED 

OUTPUT 

PHI CRIT: 

1.000 

PHI . 

CREATED 

OUTPUT 

BETA CRIT: 

1.000 

BETA 

****************************** 


ELAPSED 

TIME (SEC) 

: 0.10 




5.000 PHI 

3.200 BETA 
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******* ****** 
****** 

****** ^design gain matrix for plant ****** 

****** 

****** *********** 

DESIGN PARAMETERS : 


ALPHA = 1 . 00000E-03 

CONTROL VARIABLE C 
AIL 
RUD 

CRITERIA VARIABLE 
PHI_CRIT 
BETA CRIT 


CONTROL WEIGHT (R) 
1.0000 
2.0000 

CRITERIA WEIGHT (Q) 

4.0000 

1.0000 


STEADY STATE RICCATI SOLUTION 


1 15.95 

2 - 2.100 

3 -6.876 

4 -5.891 

5 9.802 

SAMPLING TIME : DELTA 


2 3 


-2.100 -6.876 

1.105 4.050 

4.050 15.73 

0.7099 2.397 

-0.3128 -2.262 

= 0.0000 


4 5 

-5.891 9.802 

0.7099 -0.3128 

2.397 -2.262 

3.222 -3.365 

-3.365 14.39 


******** 

******** 


EIGENVALUES OF A + B*G 
OPTIMAL CL EIGENVALUES 




COUNT 

1 

2 

3 

4 

5 


REAL PART IMAG PART 

-0.9492 0.0000 

-1.574 1.453 

-1.574 -1.453 

-4.442 2 . 532 

-4.442 -2.532 


DAMPING FREQ (RAD/S) 


1.000 0.9492 
0.7346 2.142 
0.7346 2.142 
0.8688 5.113 
0.8688 5.113 


FREQ (HZ) 

0.1511 
0.34 
0.34 
0.81 
0.81 


FEEDBACK GAIN MATRIX 


BETA 

AIL 5.470 

RUD -3.065 

ELAPSED TIME (SEC) : 


P 


PHI 


-2.605 -9.524 

-0.5127 -2.059 

1.17 


r BETA_INT 

-2.000 1.093 

2.056 -2.722 
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c********* 


****** 

****** 

****** 


^CREATE CONTROLLER FOR PLANT 


****** 

****** 

****** 



***** FULL state feedback controller ***** 
***** N0 model following ***** 

***** DIRECT f.b. states to plant ***** 

BETA P PHI R 

k * * 3 


BETA INT 


ELAPSED TIME (SEC): 0.50 



TASK 6 

k * * * 

****** ****** 

****** *FORM PLANT + CONTROLLER ****** 
****** ****** 



***** FULL STATE FEEDBACK CONTROLLER ***** 
***** MO MODEL FOLLOWING ***** 

k * * * * j 


ELAPSED TIME (SEC): 2.16 
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n***r*********************** ****** 

****** *PRINT CLOSED-LOOP SYSTEM ****** 
****** *************** 






***** CLOSED LOOP SYSTEM *** **^ ^**** ******************** *** 


SAMPLING TIME : DELTA 


0.0000 




BETA 

P 

PHI 

R 

BETA INT 


BETA P 


PHI 


BETA INT 


-0.2401 
-4.579 
0 . 0000 
7.526 
1.000 


5 . 9 64 8E-02 -4.8458E-02 -0.9132 -0.1075 

_25 99 2.8649E-02 -3.112 

0.0000 8 . 7813E-02 0.0000 

2.215 -3.809 4.611 

0.0000 0.0000 0.0000 


-8.932 
1.000 
0 . 4025 
0.0000 


===== B 



AIL 

RUD 


BETA 

2 . 5320E-03 

4 . 0504E-02 


P 

2.284 

2.060 


PHI 

0 .0000 

0.0000 


R 

0.1228 

-1.645 


BETA_INT 

0 . 0000 

0.0000 


===== c " 

BETA 

P 

PHI 

PHI CRIT 

0.0000 

1.000 

5.000 

BETA CRIT 

2 . 960 

5 . 9648E-02 

-4 . 8458E- 

AIL+ 

5 . 470 

-2 . 605 

-9.524 

RUD + 

-3.065 

-0.5127 

-2.059 

AIL= 

5 . 470 

-2 . 605 

-9.524 

RUD= 

-3 .065 

-0.5127 

-2.059 


BETA INT 


8 . 7813E-02 
-0.9132 
- 2.000 
2.056 
- 2.000 
2.056 


0 . 0000 
3 . 892 
1.093 
-2 .722 
1.093 
-2 .722 


===== d 


AIL 


RUD 


PHI CRIT 

0 . 0000 

0.0000 

BETA CRIT 

2 . 5320E-03 

4 . 050 4E 

AILE 

0.0000 

0.0000 

RUD+ 

0.0000 

0 .0000 

AIL= 

1.000 

0.0000 

RUD= 

0 . 0000 

1 . 000 

ELAPSED TIME 

(SEC) : 

0 . 12 
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TASK 8 

******* 

‘COMPUTE EIGENVALUES OF CLOSED-LOOP SYSTEM ****** 

****** 
************ 



SAMPLING TIME : DELTA = 0.0000 


******** EIGENVALUES OF CLOSED-LOOP SYSTEM ********** 


COUNT REAL PART IMAG PART 


1 -0.9492 0.0000 

2 -1.574 1.453 

3 -1.574 -1.453 

4 -4.442 2 . 532 

5 -4.442 -2.532 

ELAPSED TIME (SEC): 0.18 


DAMPING FREQ (RAD/S) 


1.000 0.9492 

0.7346 2.142 

0.7346 2.142 

0.8688 5.113 

0.8688 5.113 


FREQ (HZ) 

0.1511 

0.3409 

0.3409 

0.8138 

0.8138 



****** 


****** 


****** 

****** 


*MPAC WRITE CLOSED-LOOP SYSTEM ****** 

****** 



*** MODEL WRITTEN TO FILE: CLOSED_LOOP . MDL *** 
ELAPSED TIME (SEC): 0.32 

TOTAL JOB ELAPSED TIME (SEC): 10.24 
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SIPE TOOLBOX ARCHITECTURE 



KNOWLEDGE DATA MANAGEMENT / TYPES 

BASE 

RULES POLYNOMIAL NUMERIC 




